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We study the emission of quasi-particles in the scaling limit of the ID Quantum Ising chain at the 
critical point perturbed by a time dependent local transverse field. We compute exactly and for a 
generic time dependence the average value of the transverse magnetization, its correlation functions, 
as well as the statistic of both the inclusive and exclusive work. We show that, except for a cyclic 
perturbation, the probability distribution of the work at low energies is a power law whose exponent 
is universal, i.e. does not depend on the specific time dependent protocol, but only on the final value 
attained by the perturbation. 
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The study of out-of-equilibrium dynamics of many- 
body quantum systems has gained renewed interest over 
the past decade. This is mainly a result of impressive 
experimental advances in the field of cold atoms pi , ex- 
emplified by the observation of the collapse and revival 
of a system driven across the Mott-superfluid transition 
[2]. In many cases a common way of taking a system 
out of equilibrium is the so called quantum quench, i.e. 
a variation in time of a parameter g of the Hamiltonian 
describing the system. The resulting dynamical response 
can be probed in different ways, i.e. by looking at the 
evolution of correlation functions [3] or entanglement en- 
tropies [31 |S], or adopting a thermodynamical point of 
view and studying the statistic of the work done on the 
system [7] , the change in entropy [B] or the energy distri- 
bution after the quench [3. 

Quantum quenches can be either global or local in 
space. In both cases a change of the system parameters 
causes the emission of quasiparticles carrying correlations 
that travel across the system with a certain velocity v, 
whose maximum value for a quantum many-body system 
is given by the Lieb- Robinson bound [10] . The difference 
between the two cases lies in the fact that while in the 
global scenario the emission happens everywhere, in the 
local case it is restricted to the point where the quench is 
performed, behaving as a sort of "quantum antenna"[n]. 
Qualitatively, the result of this mechanism is the appear- 
ance of the so-called "light-cone" effect ^ , which has been 
very recently experimentally observed in bosonic ultra- 
cold gases [TT]. In the case of a local quench, on which 
we will focus in the following, the ligh-cone effect implies 
that the time lapse to see the effect of a quench at a dis- 
tance X from the antenna is t ~ x/v, the time it takes 
for excitations to reach that point [H]. The effects of a 
quench are particularly strong at a critical point where 
the excitation spectrum is gapless. In this case, a con- 
formal field theory (CFT) describes well a large class of 
Id quantum systems for sufficiently large distances and 
time scales (see for example [l3]). 

Independently on the protocol considered, the semi- 



classical picture above suggests that a local quench will 
in general produce a correlation front propagating at the 
velocity v of the quasiparticles emitted, without refer- 
ence to either the energy distribution of the excitations 
created or to the profile of the average values of physical 
observables. These quantities are however important to 
characterize how much energy is transferred to the sys- 
tem by the local perturbation and what is the form of the 
signal that propagates through the sample. These issues, 
which are in turn relevant for applications of time depen- 
dent quantum protocols in quantum information [9, 2 J 
and quantum optimization problems [14] , were hardly ad- 
dressed in the literature where most of the local quenches 
considered are "abrupt"; i.e., the parameter is suddenly 
changed from an initial value go to a final one gj |15| . 

In this Letter we consider a generic time dependent lo- 
cal quench of the transverse field in a quantum Ising chain 
at the critical point. Describing the system in its scaling 
limit, we solve exactly the dynamical problem, comput- 
ing the average transverse magnetization, its correlations 
and the statistics of the work done for a generic protocol. 
In doing so, we show that the behavior of the probabil- 
ity distribution function of the work at low energies is a 
power law whose exponent does not depend on the de- 
tailed time dependence of the protocol considered, but 
just on the final value of the local transverse field (with 
the exception of the case in which this value is zero) . We 
finally comment on the extension of these results to more 
general systems. 

Before entering the technical details of our analysis 
we will summarize and discuss the main results ob- 
tained. We will consider a quantum Ising chain, H — 
— 1/2^- crfcrf^.;^ + 5^fi where crf'^ are the longitudi- 
nal and transverse spin operators at site i and g is the 
strength of the transverse field. Setting the global trans- 
verse field g = 1, the system is at its quantum criti- 
cal point: here the mass m = 1 — g oi quasiparticle 
excitations vanishes and quasiparticles propagate freely 
through the system at the "speed of light" v ~ 1. In the 
following, we will be interested in studying the physics of 



a local change of the transverse field away from g = 1. In 
order to analyze this problem, we describe the system in 
the scaling limit by its corresponding CFT [16], locally 
perturbed by a mass term 



P,{w) 



H, 



t- -- j dx [ipd^ip - ipd^ip] + im{t)ipip\^^g, (1) 



where ip and ip are two Majorana fermionic operators, 
so that {Lp{x),Lp{x')} = {(p{x),if{x')} — S{x — x'), and 
m{t) — for all t < U), with tg arbitrary initial time, in 
such a way as to have the system in its ground state until 
t = to- Let us start by characterizing, for a generic time- 
dependent protocol, the energy distribution of quasipar- 
ticles emitted. In order to do so, imagine performing 
two measurements of energy, one before and one after 
the quench. The resulting energy difference is the so- 
called inclusive work w done on the system [171, which, 
for a system out of equilibrium, is a stochastic quan- 
tity characterized by a probability distribution Pi{w) 
[TTHTn]. We have computed its characteristic function 
Gi{u) = J dw e^^'^ Pi{w) for a generic quench protocol, 
obtaining 



Giiu) = exp 



^" J —oo J —oo 

a-i{t-t') 



dt'dtra{t)dt'm{t') 



lo. 



a — i{t — t' + u) 



(2) 



where r is the time at which the final energy measure- 
ment is performed and a is the ultraviolet cut-off of the 
theory. The work done on the system depends on the 
derivative with respect to time of the protocol m{t) cho- 
sen, matching with the expectation that this quantity 
should be related to the speed at which the quench is 
performed. 

Let us now show that the form of Pi{w) for small w, 
giving us information of the energy transmitted to the 
system, is independent on the specifics of the protocol 
employed, i.e. universal. For this sake, we have to ana- 
lyze the asymptotics of Qi for large u (much larger than 
the time scale of the protocol) . When TO(r) ^^ we get 
Qi{u) ^ (— iu)~™('^' /^^ , corresponding to 



P^(W) 



(3) 



As anticipated, the distribution Pi{w) displays an edge 
singularity at small w with an exponent independent on 
the details of the protocol m{t) chosen but just on its 
amplitude. In particular for small quenches (rn{T) < 27r) 
there is a power law divergence, while for large quenches 
Pi{w) vanishes with a cusp. This matches with the natu- 
ral physical expectation that the more the parameter in 
the Hamiltonian is changed, the smaller the probability 
of doing very small work on the system is. 

On one hand the independence of the low energy be- 
havior of the distribution of the work from the details 
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FIG. 1: (Color online) (a) Probability distributions Pi{w) for a 
nonmonotonic protocol (Blue, Full), i.e., a series of sudden 

quenches and a sudden quench (Red, Dashed) ending at the same 
value of m and shown in the inset, (b) Logarithmic plot of 
Pi(l/w) for the same protocols as before. We take q = 1. 



of the protocol can be a natural expectation in the case 
of monotonic protocols, since when u is large they all 
look like sudden quenches, making the details of how the 
final value of m is reached irrelevant. The result we ob- 
tained is, however, more general: it holds independently 
on the shape of the protocol, no matter what happens 
before the end, and therefore even in cases of nonmono- 
tonic protocols, where the former similarity is not true 
any more. We also note that, in contrast with the case 
of global quenches, where in the thermodynamic limit 
the spectral weight of the distribution is concentrated in 
a peak at high energies, the low energy part of the dis- 
tribution of the work can retain a considerable spectral 
weigth. This means that the power law behavior can be 
observable. The example of Fig. [l] clarifies both the is- 
sue of nonmonotonicity and observability. 



In Fig. 1 1(a) 

Pi{'w) is shown for a non monotonic protocol and a sud- 
den quench to the same final value of m (see the inset). 
One can see that in both cases the low energy part has 
a considerable spectral weight. From Fig. 1(b) instead 
one can see that the two protocols at low energy indeed 
behave as a power law with the same exponent. 

Cyclic protocols with 77i(t) = are an exception to the 
scenario above, since in this case the asymptotic behavior 



becomes Qi{u) ~ exp 



;fe' 



where k and k' are two 



constants depending on the specific form of m(t). In this 
case Pi{w) will have a delta- function peak, S{w) with 
a nonuniversal amplitude, plus a regular part vanishing 
linearly. This means that for each cyclic protocol there 
is a nonzero probability of not doing work on the system. 
The absence of the delta function for finite m{T) is in turn 
a consequence of the Anderson orthogonality catastrophe 

m- 

The possibility to infer the behavior of physical quan- 
tities from some gross features of the protocol m{t) is 
also observed for the transverse magnetization and its 
correlation function. The result for the average value of 
the transverse magnetization A4{x,t), represented in the 



scaling limit by the operator 2iip{x,t)ip(x,t), is 



{M{x,t)) = - 



2|x| 
7r(4a;2 + a^) 



sin(r7i(i- |x|)) . (4) 



We clearly observe that the protocol m{t) leads to the 
propagation at the velocity of light (which has been taken 
equal to 1) of two identical magnetization signals, one to 
the left and one to the right of the origin. The strength of 
the signals ,in turn, decreases with the distance as 1/x. 
The qualitative features of the traveling profile can be 
easily extracted. For example in Fig. 2(a) we have taken 



three zeros, since m{t) crosses the values n, 2it and Stt, 
and will have a positive tail, since it asymptotically ends 
at a value between Stt and Att. One can check from Fig. 



m{t) = 10(1 — e *)8(t). The traveling profile will have 



2(b) that these are indeed the features of the generated 
profile. This simple understanding can be used to design 
protocols m(t) producing a profile with the desired fea- 
tures. For instance we show in Fig. |2(c) a protocol that 
produces six positive wave-packets with the same width. 

Finally, let us show how to extend this analysis to more 
complex physical quantities, e.g. the connected correla- 
tion function of the transverse magnetization. At equal 
times this is given by 



{Mix, t)Mix', t))c = cos {m{t - \x\)) cos {m{t - \x'\)) 

^x\\x'\ 



sin {m{t — \x\)) sin {m{t — |a;'|)) 



a 



27r2 [(x - a;')2 + : 



27r2(4.T2 + a2)(4a;'2 -|- a^) 



2@{xx')xx' 



27r2(4a;2 + q;2)(4x'2 + a"^) n^ [{x + x')^ + a^] [{x - x')^ + a^ 



(5) 



Since we expect quasiparticle to be emitted symmet- 
rically, the correlation between two opposite points 
X and —X is of particular interest. In particular, 
the excess correlation C{x,t) = {M{x,t)M{—x,t))c — 
{M{x,to)M{-x,to))c, is given by 



C{x,t) 



1 

27r2(4a;2 



(cos(2TO(i- |x|))-l). (6) 




(a) 



(b) 



As before, one may easily design the protocol m{t) to give 
rise to a certain correlation profile. In particular, the ze- 
roes of C(x, t) are the same as those of the magnetization, 
as one can see from Figs. |2(b)| and 2(d) for specific proto- 
cols. More specifically, for every protocol m{t) the excess 
correlations are always negative and travel through the 
systems at the same speed as the magnetization, decreas- 
ing with the distance from the origin as l/x^. 

Let us now discuss the exact solution of this dynamical 
problem as well as the method used to obtain the results 
given above. In doing so we will also discuss in more 
detail the work done on the system, extending the re- 
sults of Eq. pi) to the exclusive work. Let us first sketch 
the main steps of the analysis. For every physical quan- 
tity at hand, it will be consistently useful to compute it 
by duplicating our theory [3T], i.e., introducing an addi- 
tional couple of Majorana fermions Xj X ^Iso described 
by the Hamiltonian of Eq. (IT]). Two pairs of Majorana 
fermions can be indeed used to form two Dirac fermions 
-ijjji = e-'^'/^iif + ix)/V2 and Vl = e*'^/4(^ + ix)V2, 
which can be combined trough a nonlocal transforma- 
tion [32], defining i'+ix) = {iPr[x) -\- V'l(— a;))/-\/2 and 
tp-{x) = {ipii{x) — 'ipL{—x))/V2i, which are described by 




(c) 
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FIG. 2: (Color online) Examples of magnetization and 

correlations profiles (right) for some specific protocol m{t) (left). 

t = 10 in (b) and t = 30 in (d). a = in all the cases. 



the Hamiltonian 



m 



i dx 



'4'Lda;ip- - il}\dx%l)+ 



-m{t) 



ip^i/j^ — '4r_'4!^ 



(7) 

This Hamiltonian describes two decoupled chiral 
modes completely characterized by the single parti- 
cle Hamiltonians -ff+._ = ^idx ± 5{x)m{t). The 
latter gives for the fields the equations of motion 



U=o 



[idt±idx]ip+-{x,t) = ±S{x)m{t)i:+-{x,t), with the 
initial condition that •ip+^-{x,to) are free fermionic op- 
erators. These equations describe the scattering of chiral 
modes on a time-dependent delta potential and their so- 
lutions can be written as 

V'+,-(x,i) = eT™(*-l-l)^(±-)V+,-(a:Ti,to)- (8) 

The computation of the average value of the magneti- 
zation and of its correlation functions proceeds now by 
expressing the operator 2i (pip — i(pip+ixx in terms of the 
^+ and "0- operators, by applying the transformations 
given above, and then by using the mode expansion of 
the fermionic operators, taking into account that the av- 
erage is taken on the Dirac sea of both fermions + and — 
(all modes with negative momentum occupied, all states 
with positive momentum free). The same procedure can 
be applied to compute the correlations. 

The probability distribution of the work can in turn 
be computed by considering that for each realization of 
an out-of-equilibrium protocol the work w is given as a 
difference of the outcomes of two measures of the en- 
ergy, at the initial time t — to and at the final time 
t — T [17) . The final energy is measured with respect 
to the full Hamiltonian Ht, or to the initial Hamil- 
tonian Ht„, excluding the forcing term, depending on 
whether one is interested in the inclusive or exclusive 
work done. For example, in the inclusive case we have 

PiM = J2n,7n^i'^- i^nir)- Ejn{to)))p{n\m,T)p,n, 

withp(n|m,r) = |(^„(r)|[/(T, to)|V'm(io))P, and p™ = 
|(V'm(io)l^(^o))l; where |$(io)) is the initial state of the 
system, f/(t,io) is the evolution operator from tg to t, 
and the equation Ht\ipi{t)) — Ei{t)\'ipi{t)) is valid. In the 
exclusive case the definition of Pe{w) is the same if En{T) 
and ipniT) are replaced by i?„(to) and ipn{to)- 

As shown in Ref. ^7\ the characteristic functions 
Gi,e{u) = J dw e™'^Pi^e(w') Contain full information about 
the statistics of w and can be written as a two-time cor- 
relations function, 



5i,e(w) 



iuHi 



'oe" 



-iuHtg\ 



(9) 



where H^f 



W {T,to)HT^toU{T,to) is the Hamiltonian 
used in the second measurement in the Hcisenberg rep- 
resentation at time r. 

To compute such quantities we use again the trick due 
to Itzykson and Zuber j21| which consists in considering 
Gi,e{u)'^- In order to do so we bosonize the Heisenberg 
representation at time r of the Hamiltonian (Is]) with the 
usual formula ^±(2;,*) == 1/V2^e=^'^^'^±(^'*\ getting 



iff = / dx {dx(t>+{x,T)f + (9,0_(x,t))V 



{dx(j)+ - dx(f>-)\^^g^ ■ 



m{T) 



(10) 



Then we observe that, apart from an irrelevant constant, 
we can write H^ = UlHtJJi, or for the sake of computing 



the exclusive work, H^ = UlHt„Ue, with Ui = e'/^^, 
Ue = e*/v^^, where 



A- 



and 



dy {(j)+{v - T, 0) + 0_(r - y, 0)) dym{T - y), 

(11) 



i? = i + m(r)(<^+(-r,O) + 0„(T,O)). (12) 

Then from Eq. ^ we get 

Gl = (^ee^^'oW.ee-™^*'') ^ {ulu^Au)), (13) 

where the average is taken on the bosonic vacuum and 
lAi^u) means that the bosonic operators are evolved at 



time u with the free bosonic Hamiltonian Ht^. Eq. (13) 
can then be computed with standard methods. 

While the result for the statistics of the inclusive work 
has been previously anticipated, for the exclusive work 
we obtain 



Ge{u) = Gi{u) X CXp 



m{T) I dtm{t) log 



—^ ( m(T)2 log —^— 
47r^ \ a ~ tu 

a' + JT- tf 

[a — iuY + (''' — t)"^ 



(14) 



When the final value of the mass m{T) is zero the inclu- 
sive (Eq. [2]) and exclusive (Eq. 14) characteristic func- 
tions coincide. Let us now briefly analyze the asymptotic 
behavior oiQ^^iu) when mir) 7^ 0. One obtains in general 



exp 



fci 



;^ 



, with fci and k'l depending on m(i), 

implying that Pe{w) consists of a delta function at zero 
as well as a regular part linearly vanishing as it; — > 0. 

Finally we briefly discuss the possibility of extending 
the result we found for Pi{w) at low energy. Given that 
a power law should be expected from the orthogonality 
catastrophe, if we make an analytic continuation u — >■ ii? 
in Eq. ([9]), we can interpret the characteristic function as 
a partition function of the classical correspondent model 
on a strip of thickness R, with a line defect stretching 
between the two boundary states U{T,to) |$(io))- The 
behavior at large R, that will determine the behavior at 
small w of Pi{w), is expected to depend on the renor- 
malization group flow of the boundary states as well as 
of the final Hamiltonian. For a local quench, we expect 
the boundary state to flow back to the critical one, since 
a non-extensive number of defects has been generated by 
the protocol. On the other hand, the defect could be rele- 
vant, marginal or irrelevant. In the case of a marginal de- 
fect (as the one explicitly considered above) the exponent 
should depend on the final strength of the defect, while 
for a relevant perturbation we do expect the exponent to 
be completely independent on the quench performed and 
equal to c/8 — 1, where c is the central charge, coming 
from the effect of a line of defect in a generic GET [21] . A 



confirmation of tfiis last statement can be observed in the 
case of sudden quenches (i.e., the Fermi edge problem) 
in Luttinger liquids [23] . 

In conclusion, we characterized the signal propagat- 
ing through the system and the energy transmitted in 
a generic time-dependent local quench of the transverse 
field in a quantum critical system, the quantum Ising 
chain at the critical point. By solving exactly the prob- 
lem in the scaling limit, we have shown that the work 
done on the system at low energies does not depend on 
the details of the protocol, but just on its amplitude. The 
independence of the physics on the fine details of the pro- 
tocol ■m{t) is observed also in the traveling signal of the 
magnetization profile and on its correlations. 
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